The detailed mechanism of the 77 production in pp scattering up to the 
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Contrary to very early beliefs, the experimental cross section data for the r) production in proton- 
proton scattering are well described if n and only r\ meson exchange diagrams are used to calculate 
the Born term. The inclusion of initial and final state interactions is done in the factorization 
approximation by using the inverse square of the Jost function. The two body Jost functions are 
obtained from the S matrices in the low energy effective range approximation. The danger of double 
counting in the pr/ final state interaction is discussed. It is shown that higher partial waves in meson- 
nucleon amplitudes do not contribute significantly bellow excess energy of Q = 100 MeV. Known 
difficulties of reducing the multi resonance model to a single resonance one are illustrated. 

PACS numbers: 13.60.Le, 13.75.Cs, 14.40.Ag, 25.40.Ep 



I. INTRODUCTION 



For a few decades the only way of detecting the rj 
meson signal in the intermediate energy proton-proton 
scattering was to perform the pionic three-prong experi- 
ments and finding the characteristic invariant mass. In 
the CERN/HERA report Q a compilation of all mu- 
tually normalized and otherwise adjusted experimental 
data for the pp — > pprj process from Ti a b = 2 GeV to 
Tiab = 13 GeV have been presented. Yet, first theoreti- 
cal models appeared at the beginning of the 90s, 
accompanied by the first published experimental results 
p| . The common denominator of these first models is a 
similar reaction mechanism: the r\ meson in two proton 
collisions is produced when an intermediate meson, 
emitted in the meson-production vertex, interacts with 
a proton forming an N* resonance which decays into r\ 
meson and proton. Everything else varies from model 
to model: the number and type of intermediate mesons 
(tt, r\, p, to, cr, ...), the type of excited state (a resonance 
formation or some re-scattering in addition), etc. The 
contribution of p meson exchange is recognized as a 
dominant one, or at least of equal importance as the 
7r meson exchange. The influence of r\ exchange has 
been considered as insignificant. The higher order terms 
expressed through different forms of initial and final 
state interactions had to be included, and it has been 
shown that they play a crucial role for the reliability of 
the calculation. The agreement of these models with the 
experimental data has been obtained by adjusting the 
free coupling constants. The new measurements have 
soon appeared 0, 0, H,EL an d the old models have been 
improved accordingly 0, Ull IT2I Il3l ] . The number and 
type of exchanged mesons still varied. 

A model has been proposed [l4| where the meson- 
nucleon partial-wave T-matrices obtained in the 



multi-resonance, coupled channel and unitary model 
|15| are used in the 77-meson production vertex instead 
of individual N* resonances. The discussion which 
intermediate mesons dominate the process is reduced to 
calculating the pp— > pprj process with the predetermined 
relative strength of the two-body eta-meson production 
T-matrices for various opened channels. Contrary to 
former statements, it has been shown that a good agree- 
ment with the experimental results for the total cross 
section is achieved using only tt and rj contributions, 
without any need to invoke heavier mesons. 

New low energy experiments reporting angular cor- 
relation measurements have been done at Jiilich [l^ . Il7| , 
so theoretical focus has moved from evaluating total 
cross section to more detailed calculations |l8j. 

This article offers an improvement of the previous 
model 14] and the extension of the examined energy 
range. Basic conclusions of the first model are confirmed: 
the 7T and rj meson exchange contributions are sufficient 
for the Born term so there is no demonstrable necessity 
to introduce other meson contributions; the initial and 
final state interactions are essential for the shape and 
size of the cross section. The new, formerly unreachable 
conclusions presented in this article are enabled by a 
more transparent treatment of the initial and final state 
interactions. This model reasonably well describes the 
total cross section for a wide energy range from very 
near threshold to Ti a b = 4.5 GeV. 



II. FORMALISM 

The ingredients of the model are presented in Fig. ^ 
The differential cross section of the full process is given 
as: 
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where T is a standard flux factor and d$ is the three- 
body Lorentz invariant phase space. The three-particle 
vertex V x is completely determined by the Bonn poten- 
tial parameters [19j, modified by keeping the momenta 
in the propagator and form factors relativistic. Pseu- 
doscalar meson-proton coupling constants and cut-off pa- 
rameters, used in this calculation, are identified with the 
corresponding on mass shell meson-nucleon values from 
the same reference; 9ttNN 

/Ait = 13.6, g% NN /4n = 3, 
= f .3 GeV and A„ = f .5 GeV. 
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FIG. 1: Ingredients of the model. ISI and FSI represent the 
initial and final state interactions. V x is the Bonn vertex, 
while the £ xv is the rj production amplitude. 



Eta emission vertex E xr) is given by xp — > r/p scattering 
amplitudes from our partial wave analysis |l5j . Off-shell 
amplitudes have been constructed by keeping the multi 
resonant T-matrices on shell and extrapolating partial 
wave projectors, following the recipe given in [20j . 

The initial and final state correction factors have roots 
in a threshold transition from the initial two proton 33 Pi 
state to the final two proton 31 So state with the r\p sub- 
system in the s-state. 

The final state correction factor is given in the factor- 
ization approximation using the two body Jost functions 
in the low energy effective range approximation |2lj | : 

f _ 1 1 (2) 
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where two body relative momentum is defined as pab — 
\Pa ~ Ps|/2 with A and B as pairs of the final state par- 
ticles. The coefficient n is a part of our model. Allowing 
the possibility that the scattering length a and effective 
range r are complex numbers, the two body correction 
factor used for the evaluation of final state contribution 
- Eq. (J2J is, in the s-wave low energy effective range 
approximation, generically given as: 

1 = P 2 -2plma + + \a + \ 2 
\J{p)\ 2 p 2 + 2plma_ + |a_| 2 ' 1 ' 

where a± = ± ± y/£ + £. 

To calculate the Jost functions we have used the follow- 
ing pp scattering length and effective range parameters 



for the two protons in the 31 So state: Oq P = 7.8 fm and 
rf = 2.79 fm P. 

For the r\p S\\ subsystem, the s-wave scattering length 
is given by d%? = (0.717-M 0.265) fm For this article, 
the corresponding effective range rf^ 1 = (—1.574 — i 0.020) 
fm has been obtained by fitting the Sn wave T-matrix 
from the same reference. 

For the 33 Pi state, the Jost function has been esti- 
mated by using p-wave effective range relation j2l|. The 
used values are a^ p — —3.3 fm 3 and r pp — 4.22 fm -1 

M 

The initial state factor fy n is indistinguishable from 
unity for energies above rj production threshold, while fi- 
nal state correction factor exhibits the distinct variations 
in magnitude. 

III. ASSEMBLING THE MODEL 

We have calculated the Born term for this process us- 
ing the Sn meson-nucleon amplitudes of ref. |Tfij for 
all three channels: it, rj, and the "effective" two body 
channel. It turns out that the individual contributions 
of 7r and rj mesons are comparable in size and slightly 
dominated by the first one. The "effective" meson con- 
tribution is negligible [l4l Wa\. Therefore, from now on 
the third meson contribution (in which the p contribution 
is implicitly included) will be disregarded in this work. 

To determine the relative sign of the remaining two 
meson contributions, we have to obtain the qualitative 
and quantitative agreement of the model predicted total 
cross section with all experimental data available within 
the low and intermediate energy range. We do it in three 
steps: we calculate the Born term contribution (model 
A), add the pp final state interaction (model B), and fi- 
nally include the pr\ final state interaction. We offer two 
versions of the pr\ final state interaction: the full Jost 
function [n = 1 in Eq. (J2J - model C), and an arbitrary 
effective reduction of the final state interaction because 
of possible double counting in the rj meson production 
vertex when the full T-matrices are used instead of indi- 
vidual N* resonances (n — \ in Eq. (J2J - model D.) 

The agreement of the model predictions with the com- 
pilation of the experimental data points to the appropri- 
ate model. 



IV. RESULTS AND CONCLUSIONS 

A. Total cross section 

The comparison of model predictions with the compi- 
lation of experimental cross section data for different 7r 
vs. r\ relative signs is given in Fig. [21 and Fig. [3] The 
negative 7r vs. r\ relative sign, as shown in FigEl is ruled 
out. 

The inclusion of the pp final state interaction improves 
the agreement of the model with the experimental data 
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FIG. 2: The total pp — > ppr] cross section. Su partial 
wave in multi resonant model is only used. The negative it 
vs. n relative sign is chosen. The interference is constructive. 
Models A,B,C and D are described in the text. Q is the excess 
energy. 
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FIG. 3: The total pp — > ppr\ cross section. Su partial wave 
in multi resonant model is only used. The positive n vs. n 
relative sign is chosen. The interference is destructive. Models 
A,B,C and D are described in the text. Q is the excess energy. 



significantly, as can be seen in Fig. [21 Used model for the 
pi] FSI gives the overall agreement between the model 
and experiment (see dashed and full lines in Fig. |3J). 



Based on the results presented in Fig. [3] we are in favor 
of model D. However, as the precision of the measured 
data in the excess energy range of 10 MeV < Q < 100 
MeV is not adequate (normalization problems), we can 
only suspect that the double counting effects are hidden 
in the proposed formalism. 
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FIG. 4: The total pp — > ppr/ cross section when individual 
Su wave resonances are consecutively added to form the two 
body T-matrix. The resonance nomenclature is taken over 
from PDG. 



The consecutive addition of individual Su wave reso- 
nances in forming the partial wave T-matrix is shown in 
Fig. 0] Results presented in this figure confirm our re- 
peatedly made statement that a single resonance model 
drastically fails in calculating the total cross section. We 
find out that only the inclusion of all three resonances 
reproduces the measured cross section in the considered 
energy range. 

In Fig. [3] the higher partial-waves are included into 
the two-body multi resonant T-matrices with the inten- 
tion of obtaining better agreement with experiment at 
energies higher then Q = 100 MeV. The inclusion of 
higher partial-waves into the 7riV T-matrices start to be 
noticeable only at surprisingly high energies (above Q 
« 300 MeV). That indicates that the improvement of 
the present model should be directed to a more rigorous 
treatment of the initial and final state interaction (pp in 
particular) prior to the possible improvement of the Born 
term. 

Additional experiments in the intermediate energy re- 
gion (100 MeV < Q < 1 GeV) would be greatly appreci- 
ated. Proper ISI and FSI treatment, along with extend- 
ing the experimental data set, might increase the relia- 
bility of the resonant parameters, the Su m particular. 
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FIG. 5: The total pp — » pprj cross section. In addition to 
the S-wave, P and D waves are as well used in xN — > r]N 
T-matrices 



B. Differential cross sections 

We have already determined the ingredients of the 
model in such a way that the shape and size of the ex- 
perimental total cross section values are fairly well repro- 
duced. The angular distributions therefore come out as a 
prediction. Comparing the predictions of model D with 
the known experimental values of differential cross sec- 
tions represent an additional test of understanding the 
physical nature of the process. 

Since the final state consists of three particles, one 
can construct and measure several partially integrated 
differential cross sections. The pprj system in the final 
state is described by 12 variables and 7 of them are fixed 
by energy-momentum conservation and particle masses. 
The differential cross sections with respect to the 5 re- 
maining quantities can be measured. We have compared 
the predictions of our model with the experimental val- 
ues for the angular distribution da/dfl rp where O,, is the 
angle of the r\ momentum in the center-of-mass frame 
with respect to the beam axis, and the invariant mass 
distribution da/ds pp where s pp is the proton-proton in- 
variant mass. Measurements of these quantities have 
been recently reported by Calen at CELSIUS facil- 
ity for Q = 15 and 37 MeV , Moskal [H at COSY-11 for 
Q = 15.5 MeV , and Abdul-Bary [H Q = 16 and 41 MeV 
at COSY TOF (time of flight) detector. The Abdul-Bary 
|25| experiment as well gives the value of da/ds pp at the 
Q = 4.5 MeV. The differential cross sections with respect 
to the remaining two variables of interest (proton-77 in- 
variant mass and the angle of relative momentum of the 
proton-proton system to the beam axis) have been mea- 
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FIG. 6: The differential cross sections. Data are from Calen 
[3 (black squares), Abdul-Bary |2fil | (gray circles), and Moskal 
[l7| (black circles). The black full lines are predictions of 
model D of this article for excess energies 4.5 MeV, 16 MeV, 
and 37 MeV respectively. The black dashed line is model D 
prediction for 16 MeV normalized by the factor 1.35. The 
gray line is model D prediction for 41 MeV. 



sured, but are not evaluated in this article. The fifth 
variable, the azimuthal orientation, is of no significance 
in the unpolarized measurements. 

Before we start the comparison we feel compelled to 
comment the overall agreement of existing experiments. 
Out of three experiments, two of them (Calen and 
Moskal 0) report quantitative measurements, while the 
third one (Abdul-Bary [^H) gives results in arbitrary 
units only, and normalizes to Calen |9j . The two quanti- 
tative experiments are clearly in strong disagreement at 
nearby energies: Calen at Q = 37 MeV |j| vs. Abdul- 
Bary at Q = 41 MeV [2l|. It is interesting to mention 
that initially reported values of the COSY-11 experiment 
pCjf are very similar to the Calen 0. However, in the 
currently valid version of the measured data the calibra- 
tion have been increased by « 40 % at Q = 15.5 MeV, 
but remained within the error bar for the Q = 4.5 MeV 
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|17| . So, we are left with the situation that we can just 
compare the shape of differential cross sections paying 
less importance to the magnitude. Disregarding the ab- 
solute normalization we see the fair agreement of the 
shape of the differential cross section for all experiments 
at Q « 15 MeV. At Q rs 41 MeV some discrepancies 
remain: the Calen [9j data show the distinct structure in 
da/dfl n which is missing in ref. |25(, and the shape of 
da/dspp is notably different in both experiments. How- 
ever, the differences remain within error bars. It is there- 
fore clear that the detailed quantitative comparison of 
theory and experiment can not be done and that achiev- 
ing the quantitative agreement awaits for the experimen- 
tal accordance. 

In Fig. El we show the results of our calculations com- 
pared to the world data collection. The full lines give the 
comparison of model D predictions with experiment. The 
present model nicely reproduces the strong dynamically 
created departure of both differential cross sections from 
the symmetric, phase space induced shape. The agree- 
ment in shape and size is good for all energies, with the 
exception of Q = 15.5 MeV value of Moskal (COSY-11) 
|l7j |. However, if we normalize our model predictions by 
35 % (Fig. 0- dashed line) the data are again reproduced 
very nicely. For the da/ds pp the low s pp peaking is nicely 
followed for all three energies, but the more pronounced 
structure at s pp ~ 3.63 GeV 2 in the Abdul-Bary data 



25] for Q — 41 MeV is not understood. The same agree- 
ment in general trend is achieved for the da / dVt^ but the 
Calen structure at Q = 37 MeV is not accounted for. 



V. SUMMARY 

Within the framework of the proposed model we 
show that the single resonance model, using only 
N(1535), drastically fails to describe the experimen- 
tal data. Next Sn resonance, N(1650), has to be in- 
cluded in any model to obtain the shape of the total 
cross section. The inclusion of the third, controversial, 
N(2090) Sn resonance represents a further improvement. 
The inclusion of P-waves is essential for the energies 
100 MeV < Q < 1 GeV. All these contributions are ef- 
fectively incorporated in the proposed model. The FSI 
corrections are essential to obtain the qualitative and 
quantitative agreement with experiment. The present 
way how the FSI is treated in this article, namely using 
only the S-waves for the Jost function in the low energy 
effective range approximation, is enough to understand 
the data to the present level of experimental precision. 
The improvement of the data, the differential cross sec- 
tions in particular, will require the improvement of the 
FSI treatment as well. 



[1] E. Flamino et al, CERN-HERA Report 84-01. 
[2] J. F. Germond and C. Wilkin, Nucl. Phys. A 518 (1990) 
308. 

[3] J. M. Laget and F. Wellers, Phys. Lett. B 257 (1991) 
254. 

[4] T. Vetter et al, Phys. Lett. B 263 (1991) 153. 

[5] A. M. Bergdolt et al, Phys. Rev. D 48 (1993) R2969; F. 

Hibou et al, Phys. Lett. B 438 (1998) 41. 
[6] E. Chivassa et al, Phys. Lett. B 232 (1994) 270. 
[7] H. Calen et al., Phys. Lett. B 365 (1996) 39; H. Calen et 

al, Uppsala University preprint TSL/ISV-95-0124 . 
[8] H. Calen et al., Phys. Rev. C 58 (1998) 2667. 
[9] H. Calen et al, Phys. Lett. B 458 (1999) 190. 
[10] C. Wilkin, Phys.Rev. C 47 (1993) R938. 
[11] G. Faeldt et al., Nucl. Phys. A 604 (1996) 441. 
[12] A. B. Santra and B. K. Jain, Nucl. Phys. A 634 (1998) 
309. 

[13] E. Gedalin, A. Moalem and L. Razdolskaja, Nucl. Phys. 

A 634 (1998) 368. 
[14] M. Batinic, A. Svarc and T.-S. H. Lee, Physica Scripta 



56 (1997) 321. 

[15] M. Batinic, et al, Phys. Rev. C 51 (1995) 2310; M. Ba- 
tinic, et al, Physica Scripta 58 (1998) 15. 

[16] J. Smyrski et al, Phys. Lett. B 474 (2000) 182. 

[17] P. Moskal, et al., Phys. Rev. C 69 (2004) 025203 

[18] K. Nakayama et al, Phys. Rev. C 68 (2003) 045201. 

[19] R. Machleidt, Adv. Nucl. Phys. 19, (1989) 189. 

[20] A. Konig and P. Kroll, Nucl. Phys. A 356 (1981) 345. 

[21] M. L. Goldberger and K. M. Watson: Collision Theory, 
John Wiley & Sons, Inc. (1964); U. Schubert, Acta Uni- 
versitas Upsaliensis 5 (1995); V. Bernard, N. Kaiser, and 
U.-G. Meifiner, Eur. Phys. J. A 4 (1999) 259. 

[22] V. G. J. Stoks and Th. A. Rijken, Phys. Rev. C 59 (1999) 
3009. 

[23] J. R. Bergervoet, et al, Phys. Rev. C 38 (1988) 15. 
[24] S. Ceci and A. Svarc, Few-Body Systems Suppl. 14 (2003) 
311. 

[25] M. Abdul-Bary et al, Eur. Phys. J. A 16 (2003) 127. 
[26] P. Moskal et al, PiN Newslett. 16 (2002) 367. 



